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Abstract. We discuss the coarse-grained level density of the Henon-Heiles 
system above the barrier energy, where the system is nearly chaotic. We use 
periodic orbit theory to approximate its oscillating part semiclassically via 
J> | Gutzwiller's semiclassical trace formula (extended by uniform approximations 

^n ■ for the contributions of bifurcating orbits). Including only a few stable and 

. unstable orbits, we reproduce the quantum- mechanical density of states very 

accurately. We also present a perturbative calculation of the stabilities of two 
infinite series of orbits (R n and L m ), emanating from the shortest librating 

■ straight-line orbit (A) in a bifurcation cascade just below the barrier, which 
at the barrier have two common asymptotic Lyapunov exponents xr an d XL- 
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The two-dimensional Henon-Heiles (HH) Hamiltonian 



H HH = T + V H h(x, y) = \ {pi +p 2 y ) + \ (x 2 + y 2 ) + a (x 2 y - y 3 /3) (1) 



was introduced to describe the mean gravitational field of a stellar galaxy. It 
describes an open system in which a particle can escape over one of three barriers 
with critical energy Eb ar = l/6a 2 and has meanwhile become a textbook example 
15 °f a system with mixed dynamics reaching from integrable motion (for 
E — * 0) to nearly fully chaotic motion (for E > Eb ar ). Scaling coordinates and 
momenta with a causes the classical dynamics to depend only on the scaled 
energy e = E/Eb ar = 6a 2 E; the barrier energy then lies at e = 1. 

The Hamiltonian Q has also been used jS] to describe the nonlinear nor- 
mal modes of triatomic molecules, such as H^~, whose equilibrium configuration 
has -D3 symmetry. Although this model may no longer be quantitative for large 
energies, it can qualitatively describe the dissociation of the molecule for e > 1. 

In this paper we discuss the coarse-grained level density of the HH Hamilto- 
nian Q above the barriers, calculated both quantum-mechanically and semiclas- 
sically using periodic orbit theory. Since the potential Vhh m © goes asymp- 
totically to —00 like — r 3 (r 2 = x 2 + y 2 ) in some regions of space, the quantum 
* E-mail address: matthias.brack@physik.uni-regensburg.de 
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spectrum of (^Q) is strictly speaking continuous. However, for sufficiently small a 
there are quasi-bound states for E < Eb ar whose widths are exponentially small 
except very near E^ ar . For semiclassical calculations of the HH level density for 
e < 1, we refer to earlier papers [HI d - I n we have calculated the complex 
resonance energies E m — iF m by the standard method of complex rotation, di- 
agonalizing in a finite harmonic-oscillator basis. The level density is, after 
subtracting the non-resonant part of the continuum, given by 



Ag(E) =-IlmV 

TV ' 



1 



7T ^ E-E m + ir m /2 ' 

m 

We define the coarse-grained level density by a Gaussian convolution of ((21) over 
an energy range 7 

1 /*oo 

Ag,(E) = -±= / Ag(E') e -(^') 2 /7 2 dE >, (3) 

which can be done analytically [H]. Its oscillating part, which describes the gross- 
shell structure in the quantum-mechanical level density, is then given by 

5g qm (E) = Ag, / (E)-Ag(E), (4) 

where Ag(E) is the smooth part of © which we have extracted by a complex 
version [S] of the numerical Strutinsky averaging procedure [S] . 

Semiclassically, the quantity 5g(E) can be approximated by Gutzwiller's trace 
formula [TJU, which for a system with two degrees of freedom reads 



Sg sd (E) = ±Y Tp ° (E) e -hTME)/2H] 

nh^ rpo ^\TrM po (E)-2\ 
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Spo(E) _ 7T 

h 2 Upo 



(5) 



The sum goes over all isolated periodic orbits labeled 'po', and the other quan- 
tities in ® are the periods T po and actions S po , the Maslov indices o po and 
the repetition numbers r po of the periodic orbits. M po (E) is the stability matrix 
obtained by linearization of the equations of motion along each periodic orbit. 
The Gaussian factor in (JSJ) is the result of a convolution analogous to (0) ; it sup- 
presses the orbits with long periods and hence yields the gross-shell structure in 
terms of the shortest periodic orbits, hereby eliminating the convergence prob- 
lem characteristic of non-integrable systems This use of the trace formula 
to describe gross-shell quantum effects semiclassically has found many applica- 
tions in different fields of physics (including interacting fermion systems in the 
mean- field approximation; see [I] for examples). 

The shortest periodic orbits of the classical HH system Q have already been 
extensively studied in earlier papers ^2 El D2j In [E] we have calculated all 
relevant orbits and their properties from the classical equations of motion and 
computed the quantity 5g sc i (E) in ® . Some of the shortest orbits are shown in 
the left part of Fig. ^ all evaluated at e = 1 (except for T2 which is evaluated at 
e = 1.1). Note that due to the D% symmetry of the HH potential, the orbits A5, 
B4, T2 and L/6 (as well as all orbits R n and L m bifurcating from A, see below) 
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Figure 1. Left: Contours of the HH potential and some of its shortest periodic orbits in the 
{x,y) plane (see text). Right: Trace of the stability matrix of the A orbit and the three pairs of 
orbits (R,5,L,6), (R7,Lg), (R,9,Lio) bifurcated from it, forming the beginnings of the 'HH fans'. 
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Figure 2. Comparison of quantum-mechanical (solid line) and semiclassical (dashed line) level 
density Sg(E) of the HH potential versus scaled energy e, coarse grained with Gaussian smooth- 
ing range 7 = 0.25. Only 18 periodic orbits contribute to the semiclassical result [Sj. 



have two symmetry partners obtained by rotations about ±27r/3. The orbit C3 
and all triplets of R m orbits have a time reversed partner each. 1 

In Fig.Elwe show a comparison of semiclassical © with quantum-mechanical 
@ results, both coarse-grained with 7 = 0.25 (units such that h = 1). At this 
resolution of the gross-shell structure, only 18 periodic orbits contribute to the 
semiclassical result; for the period-two orbit D7 which is stable up to e ~ 1.29 and 
involves to further orbits (Eg, G7) in a codimension-two bifurcation scenario, we 
have used the appropriate uniform approximation to avoid the divergence 
of the trace formula © (see [B] for details). We note that the agreement of 
semiclassics with quantum mechanics is excellent. Only near e ~ 1 there is a 
slight discrepancy which is mainly due to some uncertainties in the numerical 
extraction of Ag(E). We can conclude that also in the continuum region above 

1 We use here the nomenclature introduced in |12II13| . where the Maslov indices a po appear as 
subscripts of the symbols (B4, R5, L6, etc.) of the orbits. 
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a threshold, the semiclassical description of quantum shell effects in the level 
density of a classically chaotic system works quantitatively. 

In view of the importance of the level density close to the critical barrier 
energy e = 1 for the threshold behaviour of a reaction described by the HH model 
potential, we focus now on a particular set of periodic orbits existing at e = 1. 
The straight-line librating orbit A reaches this energy with an infinite period after 
undergoing an infinite cascade of bifurcations for e — > 1. At these bifurcations, 
two alternating infinite sequences of rotational orbits R n (n = 5,7,9,...) and 
librating orbits L m (m = 6, 8, 10, . . .) are born; their bifurcation energies e n and 
e m form two geometric progressions converging to e = 1 with a 'Feigenbaum 
constant' 5 = exp(27r/v / 3) = 37.622367...; the shapes of these new orbits are 
self-similar when scaled with y/d in both x and y direction (see ^2] for details). 
The stability traces of the first three pairs (R n , L m ) are shown in the right part 
of Fig. n As seen there, the curves TrM(e) of these orbits are nearly linear (at 
least up to e ~ 1.02) and intersect at e = 1 approximately at the same values for 
each type (R or L). For large n and m, these values were found numerically J2j 

to be TrM Lm (e = 1) ~ 8.183 and TrM Rn (e = 1) 4.183. This means that at 

e = 1, all L orbits have asymptotically the same Lyapunov exponent xl — 2.087, 
and all R orbits have the same Lyapunov exponent xr — 1-368. Based upon 
these numerical findings, we postulate the following asymptotic behaviour: 



TrM RnjLm (e) ~ 2 T 6.183 — — for e -► 1 , n, m - oo . (6) 



Here e* are the respective bifurcation energies of the orbits (e n or e m ), and the 
minus or plus sign is to be associated with the R or L orbits, respectively. The 
curves TrM Rni L m (e) thus form two 'fans' spreading out from the values 8.183 and 
—4.183 at e = 1, the first three members of each being shown in the right part of 
Fig. ^ In the following we will sketch briefly how the qualitative features in Q 
of these 'HH fans' can be obtained analytically from semiclassical perturbation 
theory. Details will be given in a forthcoming publication |15j . 

The idea is to start from the following 'separable HH' (SHH) Hamiltonian 



and to include the term ax 2 y in first-order perturbation theory. Formally, we 
multiply it by a small positive number e and write Hhh = Hq + e H\ with 
Hi = ax 2 y = u 2 v/a 2 , where u = ax and v = ay are the scaled coordinates. The 
Hamiltonian (JJJ) is integrable; an analytical trace formula for it has been given in 
0. There is only one saddle at (x,y) = (0, 1) with energy e = 1 and one librating 
A orbit along the y axis which undergoes an infinite cascade of bifurcations for 
e — * 1. From it, an infinite sequence of rational tori bifurcates, where I is 
their repetition number and k = 2, 3, . . . counts the bifurcations (and the tori). 
The v motion of the primitive A orbit (I = 1), having u^(t) = 0, is given by 




#0 = H S hh = g (Pi +P 2 y) + 2 ^ + y2) " 3 y3 



(7) 
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in terms of the Jacobi elliptic function ^H] sn(z,q) with modulus q. In (|8j). 
vi < V2 < i>3 are the turning points of the motion along the v axis, defined by 
Vhh(u=0,Vi) = e (i = 1,2,3). The tori bifurcating at the energies have the 
same v motion as the A orbit: i>t(£) = i>a(£). Their u motion is given by 



u T (t) = y/( e -e lk )/3 sin(t + 0), e > e, fc , </>G[0,2vr). (9) 

The angle <f> describes the members of the degenerate families of tori. 

According to semiclassical perturbation theory ^7j; the actions Sik of the 
tori are changed in first order of e by 

hSihif) = -e 6 H 1 (u T {t),v T (t))dt = / lk u^{t)v T (t)dt, (10) 

Jlk a Jo 

where T$P' = 2irk are the periods of the unperturbed tori 0. The integral in 
H1U|) takes the form 5iSik((f>) = A.\k + Bik cos(2<^>). In the asymptotic limit e — * 1, 
where g — * 1 and ~ ln[432/(l — e)], the coefficients Ai k and Bu- can be given 
analytically [T^l- Integrating the phase shift caused by 6iSik(<t>) over the angle 
(i.e., over the torus Tik) yields a modulation factor -A^ifc 

Mi k = ±- r e^ s >^ d(f> = e^ k M\B lk \/K) , (11) 

27T Jo 

to be inserted under the sum of tori in the trace formula for the unperturbed SHH 
system given in [7j. Replacing the Bessel function in (|11|) by its asymptotic form 
Jo(x) ~ y/2/irx cos(x— 7r/4) yields two terms for each torus Tik, corresponding to 
the two isolated orbits R and L into which it is broken up by the perturbation. 
Reading off their overall amplitudes ^4r,l in the perturbed trace formula and 
identifying them with their expression for isolated orbits given in (JSJ, i.e. equating 

4L = T ; ^ (12) 
vrft V|TrM R , L -2| 

uaing the unperturbed periods T$ for Tr^, we can determine the perturbative 
expression for the stability traces. For the first repetitions (1 = 1) they become 

TrM RniLm (e) ~ 2 =F 5.069 for e -> 1 , (13) 

\l-e lk J 

and thus have exactly the same functional form as in ©. Here &\ k are the 
bifurcation energies of the primitive A orbit; k = 2, 3, . . . labels the pairs of R m 
and L n orbits with m = 2k + 1 and n = 2k + 2, and the signs are to be chosen 
as in (JBJ). In (|13|) we have put e = 1 which is justified since even for this value 
the perturbations 5\Sik near e ~ 1 are sufficiently small. 

Although the perturbative result ()13|) contains a too small value of the con- 
stant 5.069 (instead of 6.183) by 18%, it explains qualitatively correctly the 
numerical features of the 'HH fans' in ©, in particular the linear intersection of 
the curves TrM Rj L(e) at e = 1 at two values lying symmetrically to TrM = +2. 
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